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Abstract. We present a Green's dyadic formulation to calculate the Casimir energy 
for a dielectric-diamagnetic cylinder with the speed of light differing on the inside and 
outside. Although the result is in general divergent, special cases are meaningful. It is 
pointed out how the self-stress on a purely dielectric cylinder vanishes through second 
order in the deviation of the permittivity from its vacuum value, in agreement with 
the result calculated from the sum of van der Waals forces. 
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1. Formulation of the Green's dyadic approach 

The electromagnetic Green's dyadic functions have been succcesfuUy used in many 
occasions (for an extensive view see |2] and references whitin) and can be applied to 
very complicated geometries. Their use happen to be critical in this calculation 
This approach helps us compute the vacuum expectation value of the fields rigorously; 
we show that the approach is both illuminating of the physics and unambiguous. 

1.1. Green's dyadic equations; formalism 

In a medium of constant electric permittivity e' and magnetic permeability fi' we insert 
an infinitely long cylinder of radius a with permittivity and permeability e and fi. The 
product of these parameters is different than that of the outside parameters. There are 
no real charges of any kind present in the problem, p = J = and since we work at a 
fixed frequency we can Fourier transform the electric and magnetic fields, 

E(r,t)=/ -^E(r,.;)e— *, B(r,t)=/ -^B(r,^)e— (1) 
and the corresponding Maxwell's equations are 

V X E = icjpH, V-D = 0, (2a) 

V X H = -iuje'E, V-B = 0. (26) 
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In order to write down the Green's dyadic equations, we introduce a polarization 
source P. The first equation in (j2^ and the second one in ()2aj) get then changed to, 

V X H = -iueE - iuP, V-D = - V- P. (3) 

The hnear relation of polarization source with the electric field defines the Green's dyadic 
as 

E(x) = J {dx')T{x,x') • P(x'). (4) 

Since the response is translationally invariant in time, we work with the Fourier 
transform of the dyadic at a given frequency u. We can then, by simple substitution, 
write the dyadic Maxwell's equations in a medium characterized by a dielectric constant 
e and a permeability yU§: 

V xT' -iujfi{uj)^ = _L-V X 1, V-$ = 0, (5a) 

-V X $ - iuje{uj)T' = 0, V • T' = 0. {5b) 

and where the unit dyadic 1 includes a three-dimensional 6 function, 1 = 15(r — r'). 
Quantum mechanically, these Green's dyadics give the one-loop vacuum expectation 
values of the product of fields at a given frequency u, 

(E(r)E(r')) = -r(r, r'), (H(r)H(r')) = --^^ x r(r, r') x V'. (6) 
1 1 w^/i^ 

Thus, from the knowledge of the classical Green's dyadics, we can calculate the vacuum 

energy or stress. 

Since the TE and TM modes do not separate, we cannot use the general waveguide 
decomposition of modes into those of TE and TM type||. However we can introduce the 
appropriate partial wave decomposition for a cylinder, in terms of cylindrical coordinates 
{r,e,zn: 

r'(r,r';^)= J2 / 7^ (V x z)/„(r; A:, 2;) 

m=-oo^-'» 27r I 



+ —V X (V X z)gm{r; k, uj)Xmk{6, z) i, (7a) 
ue J 

°° r°° dk ( 

$(r,r;u;) = ^ / —\{Vx%)g^{r]k,u;)xmk{0,z) 

- -^V X (V X z)fm{r;k,uj)xmk{0,z)], (76) 

where the cylindrical harmonics are x(^; z) = ^^^e"^^e''^^, and the dependence of fm 
etc. on r' is implicit. Notice that these are vectors in the second tensor index. Because 
of the presence of these harmonics we have 

1777/ ^ Q A fYtJxi 

V X z ^ f 0— = Al, and V x (V x z) ^ fiA;- 6 zd^ = Af, (8) 

r or or r 

§ In order to have divergenceless Green dyadics, we redefine the electric Green's dyadic in the following 
way, r'(r, r', ui) = r(r, r', ui) + j^<5(r — r') and $ is the magnetic dyadic. 

II For example as given in Ref. However, this is here impossible because the TE and TM modes do 
not separate. See Ref. 

% A slight modification of that given for a conducting cylindrical shell |H| ■ 
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in terms of the cylinder operator d^ = ^-^r-^ — ^. It is trivial to see that the divergence 
of ()7ap and (ff^ is zero, satisfying immediately two the the dyadic Maxwell's equations. 
Now, if we use the Maxwell equation (j56jl we conclude"*" 



9n 



9n 



and 



{dm k )fm — ^ fJ'fn 



(9) 



More elaborate work is needed to get a condition form the other Maxwell equation 
fjHffljl . Using the above we can write (jHajl as, 

r°° dk \ ^^{d, '-^^ 

-oo 277 



E 



M 



'^T—^fm - —{dm - k^Wgm \Xmk{6, z) 



°° dk 



J2 I ^\ iujfxAfgm + eMfm \xmk{d, z) + -V xl . 



1. 



(10) 



if we multiply the above by the expression J^"" J^^d9dzXm'k'{^^ apply 

It I-o. dedzx*m'AO. z)xmk{0, z) = 2'K5{k - k')6mm' , we find 

1 ~ i 
:r-Af{dm -k"^ + uj^fie)fm M,{dm - k^ + uj^fie)gm = 



1 r2iv fOO 1 

- / / d9dzx*mk{0,z){V X l)-S{r-r')6{9-9')6{z~z'), 
s Jo J~oo r 



where the delta functions are now made explicit. By dotting this expression with z we 
notice that z-Ai = and z-J\f = —dm and after a little manipulation we get to the 
fourth order differential equation: 



d V i (r- r' 6' z' 

'-'"m^ m'-m\i i ' ; ; ^ 



^M'*U{r - r')XmuKo , 
e r 



z') 



(12) 



(13) 



If we now dot it with (Vxz), we learn that a similar equation holds for Qm'- 

dmVmgUr; r\ 9', z') = -iu;Af'*U{r - r')Xm,{9\ z'), 

r 

where we have made the second, previously suppressed, position arguments explicit and 
the prime on the differential operator signifies action on the second primed argument*. 
To solve those equations, we separate variables in the second argument. 



fm(r,r') 



gm(r,r ^ 



M'*Fm{r, r'- k, iu) + -M'*Fm{r, r'- k, u) 

UJ 



Ar'*Gm{r, r- k, - iM'*Gm{r, r ■ k, uo 

UJ 



XmkK^'.z'), (14a) 

X*mr.{0\z)lUh) 



where we have introduced the two scalar Green's functions F^, Gm satisfying 

dmT>mFm{r,r') = '^^-S{r - r'), and dmT>mGm{r,r') = uj^-6{r - r'), (15) 
e r r 

while Fm and Gm are annihilated by the operator dm'Dm, 

dmVmF{r, r') = dmVmG{r, r') = 0. (16) 

+ The ambiguity in solving for these equations is absorbed in the definition of subsequent constants of 
integration. 

* The Bessel operator appears, 2?„i = dm + A^; = oj^e^ — fc^. 
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1.2. Green's dyadic solutions 



The Green's dyadics have now the form: 

°° r°° dk ( / d - P\ 1 



m=—co 



+-MJ\f'* 



dm k 
dk 



^ rOO riff ( 1 ip 

$(r,r»= / ^\--MM"'Gm{ry)-—J\fM"'Fm{ry) 



(176) 



In the following, we will apply these equations to a dielectric-diamagnetic cylinder of 
radius a, where the interior of the cylinder is characterized by a permittivity e and 
permeability /i, while the outside is vacuum, so e = /i = 1 there. Let us consider the 
case that the source point is outside, r' > a. If the field point is also outside, r, r' > a, 
the scalar Green's functions F^, G^, F', G' that make up the above Green's dyadics 
(we designate with primes the outside scalar Green's functions or constants) obey the 
differential equations (fT3j) and (fTB|) with e = = 1. The solutions to these equations 



aren: 



F' (r r') = — 



/\m\ 



-\m\ 



\m\ 



X'"^ 2\m\ \r- 



+ 



+ B'lHUX'r' 



Hm{>^'r) - — — J„(AV<)if„(AV>), 



while G'^ has the same form with the constants a'^ , b'J^ , A'^ , B'^ replaced by 
a'^, b'^, A'^, B'j^, respectively. The homogeneous differential equations have solutions 



F' (r r') = — 



-.IF 



n/|m| 



+ C^m(AV) 



_|_ 



A. 



IF 



+ B'iHU\'r') 



i\in\ 



HM'r), (19) 



a"^, etc. 



while in G'^ we replace a'^ 

When the source point is outside and the field point is inside, all the Green's 
functions satisfy the homogeneous equations (fT^ with e, fi 1, and then 
Fm,Gm,Fm,Gm, are of the same form as in equation (fT!I|) with the corresponding 
change of constants. In all of the above, the outside and inside forms of A are given by 



A 



12 



k"^ and A^ = uj'^^e — k"^ 



The various constants are to be determined, as far as possible, by the boundary 
conditions at r = a. The boundary conditions at the surface of the dielectric cylinder are 
the continuity of tangential components of the electric field, of the normal component 
of the electric displacement, of the normal component of the magnetic induction, and 
of the tangential components of the magnetic field (we assume that there are no surface 

tt For details see 13 and [H]. 
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charges or currents). In terms of the Green's dyadics, the conditions read 



v-eT' 



r=a+ 



r=a— 
r=a+ 



0, 



0, 



T' 



e 



r=a+ 



r=a- 
r=a+ 



0, 



0, 



z-r' 



r=a— 
■=a+ 



0, 



0. 



(20 a) 
(206) 



By imposing those boundary conditions, we find that the only constants contributing 
to the energy are: 

mkuj 



B: 



G _ 
m 








'G _ 




m 





XX'D 



Jm{Xa)H^{\'a)B:, 



F 

m' 



mkuj 



B: 



Bi 



IF 




iJ^ n Jm{X'a) 
A^2ii7„(A'a 

- 4(1 - efi)^MXa)HU\'a)B: 
XX'D 



G 
ml 



B. 



IF 



1 , , mkuo o , , , ^ 



G 
ml 



iG 



B: 



all in terms of BE, 



up' TT Jm{X'a) 

X^2AH^{X'a 



■^M and 




1 JmjXa) 
e Hm{X'a 



■B 



G 



(21a) 
(216) 

(21c) 
(2 Id) 
(21e) 

(21i) 



B, 



G 



D 



The denominators occurring here areff 

E = {1 - efif———Jl{Xa)H^{X'a) 



DD, 



(22 a) 
{22b) 
(22 c) 



D = eX'aJ'^{Xa)Hm{X'a) — XaH'^{X' a) Jm{Xa) , 
D = fiX'aJl^{Xa)Hm{X'a) — XaH'^^^X' a) Jm{Xa) . 

It is now easy to check that the terms in the Green's functions that involve powers 
of r or r' do not contribute to the electric or magnetic fields. So, even though we are 
not able to determine all the constants (notice that there is some ambiguity in these 
since they cannot be uniquely determined), it is not an issue since the energy will be 
well defined jSJIH]. These constants enter always in the same form and therefore their 
individual values are not relevant. As we might have anticipated, only the pure Bessel 
function terms contribute. It might be thought that m = is a special case, and indeed 
— > |ln^, but just as the latter is correctly interpreted as the limit as 
\m\ 0, so the coefficients in the Green's functions turn out to be just the m = limits 
for those given above, so the m = case is properly incorporated. 

ft The denominator structure appearing in S is precisely that given by Stratton 0,and is the basis for 
the calculation given, for example in Ref. \7\. It is also employed in an independent rederivation of the 
Casimir energy for a dilute dielectric cylinder [S]. 
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1.3. Stress on the cylinder 

We are now in a position to calculate the pressure on the surface of the cylinder from 
the radial-radial component of the stress tensor 

P=(T„,)(a-)-(T„,)(a+) (23) 

where Trr = | [e{E^ + El - E^.) + + - H^)]. As a result of the boundary 
conditions, the pressure on the cylindrical walls is given by the expectation value of the 
squares of field components just outside the cylinder, therefore 

e-l 



T 



— T 

J. r-i 



Ei 



(24) 



a+ 



where the expectation values are given by © in terms of the Green's functions. We 
obtain the pressure on the cylinder as 



P 



E 



X Km{y' 

,2 



m=—oo 



dCadka^i K'^{y')I^{y)I'^{y)y{k a - ( a fi) - K'^{y')I^{y) 



y'^e 



Pa' - Ca'e 



y 



m 



y' V 



Ca 



+ y 



12 



K'M)C{y)KUy')t^y'{k'a' - ('a'e 



Uy)Uy)KUy')y 



m 



""nil yi 2\ . 12 
— [k a n - C a ) +y fi 



y 



+ {(£^/i)}. 



(25) 



where we have performed the Euclidean rotation u ^ i(, X ^ in, and S is the rotated S. 
Here y = na, y' = n'a and the last bracket indicates that the expression there is similar to 
the one for the electric part by switching e and fi, showing manifest symmetry between 
the electric and magnetic parts. However, this expression is incomplete. It contains 
an unobservable "bulk" energy contribution, which the formalism would give if either 
medium, that of the interior with dielectric constant e and permeability fi, or that of 
the exterior with dielectric constant and permeability unity, fills all the space The 
corresponding stresses are computed from the free Green's functions which satisfy (|THjl . 
and have solutions 



1 



|m| 



TT 



2 m \ r 



2i 



(26) 



where < r, r' < oo. Notice that in this case, both F^^ and are zero. After the 
Euchdean rotation the bulk pressure becomes 

n 



pb = Tj:^)^a-)-Tj:^{a+) 



E 



m=— oo 



dCa dka{y^I'^{y)K'^{y) 



- {y' + m')Uy)Km{y) - y"Uy')K'M) + {y" + m')iM)KM)}- (27) 

This term must be subtracted from the pressure given in (123). Note that P'' = in the 
special case e/x = 1 as it should be. 
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2. Dilute dielectric cylinder 

We now turn to the case of a dilute dielectric medium filling the cylinder, that is, set 
fi = 1 and consider e — 1 as small. We can then expand the integrand in (j^^ and (jTfjl 
in powers of (e — 1). Because the expression in p5|l is already proportional to that 
factor, we need only expand the integrand to first order. The total pressure can then 
be written as: 

.2\ 



+ KZiy)Ciy)l - KZ{y)ll{y)\ ( i + ^ 1 + l^'^W^iy)Km{y) 



+ KlwUy)l f 1 + ^ ) ( 1 - - Kliy)C{y)l (i - ^ 



+ Kl{y)C{y)Uy) (l + ^ j ] + (28) 

Thus the total stress vanishes in leading order which is consistent with the interpretation 
of the Casimir energy as arising from the pairwise interaction of dilutely distributed 
molecules. Several methods to compute this integral are explain with great detail in jS] 
and in jHj . There it is shown that making use of the asymptotic expansion for the Bessel 
functions, we can numerically evaluate the integral 

P = ^(-0.007612 + 0.287168 + 0.024417 - 0.002371 - 0.000012 - 0.301590) 

327r"'a'^ 

= 0.000000, (29) 

and by introducing an exponential regulator e~^y in (j28|) we can unambiguously separate 
the two divergent terms 

327r2a4 1^32^3 81925 j ' ^ ^ 

The form of the divergences is exactly as expected fH^j. In particular, there is no 1/5^ 
divergence. How do we interpret these terms? It is perhaps easiest to imagine that 6 as 
given in terms of a proper-time cutoff, S = r/a, t ^ 0+. Then if we consider the energy, 
rather than the pressure, the divergent terms have the form ii^div = 63^ + 61-^^. Here 
L is the (large) length of the cylinder. Thus, the leading divergence corresponds to an 
energy term proportional to the surface of the cylinder, and it therefore appears sensible 
to absorb it into a renormalized surface energy which enters into a phenomenological 
description of the material system. The 1/r divergence is more problematic. It is 
proportional to the ratio of the length to the diameter of the cylinder, so it seems 
likely that this would be interpretable as an energy term referring to the shape of 
the body. In any case, although the structure of the divergences is universal, the 
coefficients of those divergences depend in detail upon the particular regularization 
scheme adopted. The nature of divergences in such Casimir calculations is still under 
active study j21 CHI CH CS] ■ In contrast, the term proportional to {e — l)^/a^ is unique. 
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The universality of the finite Casimir term makes it hard not to think it has some real 
significance. Thus, of course, it could not have been any other than that zero value 
given by the van der Waals calculations [3 UHl Hlj . 

3. Conclusion 

We have shown how the Green's dyadic formulation, modified for dielectric materials, 
exhibits a transparent way to calculate the Casimir energies of a dielectric-diamagnetic 
cyhnder and showed that in the dilute case, it coincides with that obtain by summing the 
van der Waals energies of the constituent molecules. However, the identity is not really 
that trivial, because both the van der Waals and the Casimir energies contain divergent 
contributions. This is particularly crucial when one is considering the self-stress of a 
single body rather than the energy of interaction of distinct bodies. It was nontrivial 
to show the analog for the case of the dielectric sphere JH] , and the calculation for the 
dielectric cylinder turned out to be extraordinarily difficult. 
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